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Abstract 

We carry out a comprehensive analysis of the nonminimal supersymmetric 
standard model (NMSSM) with baryon and lepton number violation. We 
catalogue the baryon and lepton number violating dimension four and five 
operators of the model. We then study the renormalization group evolution 
and infrared stable fixed points of the Yukawa couplings and the soft super- 
symmetry breaking trilinear couplings of this model with baryon and lepton 
number (and R-parity) violation involving the heaviest generations. We show 
analytically that in the Yukawa sector of the NMSSM there is only one in- 
frared stable fixed point. This corresponds to a non-trivial fixed point for the 
top-, bottom-quark Yukawa couplings and the B violating coupling A233, and 
a trivial one for all other couplings. All other possible fixed points are either 
unphysical or unstable in the infra-red region. We also carry out an analysis 
of the renormalization group equations for the soft supersymmetry breaking 
trilinear couplings, and determine the corresponding fixed points for these 
couplings. We then study the quasi-fixed point behaviour, both of the third 
generation Yukawa couplings and the baryon number violating coupling, and 
those of the soft supersymmetry breaking trilinear couplings. From the anal- 
ysis of the fixed point behaviour, we obtain upper and lower bounds on the 
baryon number violating coupling A233, as well as on the soft supersymmetry 
breaking trilinear couplings. Our analysis shows that the infrared fixed point 
behavior of NMSSM with baryon and lepton number violation is similar to 
that of MSSM. 
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I. INTRODUCTION 



Supersymmetry |TJ is at present the only known framework in which the Higgs sector of 
the Standard Model (SM), so crucial for its internal consistency, is natural. A much favored 
implementation of the idea of low energy supersymmetry is the Minimal Supersymmetric 
Standard Model (MSSM), obtained by simply doubling the number of states of the SM, and 
introducing a second Higgs doublet (with opposite hypercharge to the SM Higgs doublet) to 
generate masses for all the fermions and to cancel traingle gauge anomalies. The minimal 
supersymmetric version of the standard model leads to a successful prediction for the ratio of 
the gauge couplings with a gauge unification scale M G ~ 10 16 GeV. This has led to the idea 
that there may be a stage of unification beyond the SM. If so, then it becomes important 
to perform the radiative corrections in determining all the dimension < 4 terms in the 
lagrangian. This can be achieved by using the renormalization group equations in finding 
the values of parameters at the low scale, given their value at a high scale. Thus, considerable 
attention has recently been focussed on the renormalization group evolution || of the various 
dimensionless Yukawa couplings in the SM and its minimal supersymmetric extension, the 
Minimal Supersymmetric Standard Model. Using the renormalization group evolution, one 
may attempt to relate the Yukawa couplings to the gauge couplings via the Pendleton-Ross 
infra-red stable fixed point (IRSFP) for the top-quark Yukawa coupling ||, or via the quasi- 
fixed point behaviour [|J . The predictive power of the SM and its supersymmetric extensions 
may, thus, be enhanced if the renormalization group (RG) running of the parameters is 
dominated by infra-red stable fixed points (IRSFPs). Typically, these fixed points are for 
ratios like Yukawa couplings to the gauge coupling, or in the context of supersymmetric 
models, the supersymmetry breaking trilinear A-parameter to the gaugino mass, etc. These 
ratios do not attain their fixed point values at the weak scale, the range between the GUT (or 
Planck) scale and the weak scale being too small for the ratios to closely approach the fixed 
point. Nevertheless, the couplings may be determined by quasi-fixed point behaviour 0] 
where the value of the Yukawa coupling at the weak scale is independent of its value at the 
GUT scale, provided the Yukawa coupling at the GUT scale is large. For the fixed point 
or the quasi-fixed point scenarios to be successful, it is necessary that these fixed points are 
stable §-|7|. 

Since supersymmetry |IJ requires the introduction of superpartners of all known particles 
in the SM (in addition to the introduction of at least two Higgs doublets), which transform 
in an identical manner under the gauge group, there are additional Yukawa couplings in 
supersymmetric models which violate H baryon number (B) or lepton number (L). In 
MSSM, a discrete symmetry |§ called R-parity (R p ) is invoked to eliminate these B and 
L violating Yukawa couplings. However, the assumption of R p conservation at the level of 
the MSSM appears to be ad hoc, since it is not required for the internal consistency of the 
model. Therefore, the study of the renormalization group evolution of the dimensionless 
Yukawa couplings in the MSSM, including B and L (and R p ) violation, deserves serious 
consideration. Recently considerable attention has been devoted to this question [ID, IT 



It has been shown that the only stable infrared fixed point of MSSM with baryon and 
lepton number violation is the one where the top-, bottom-quark Yukawa couplings and the 
B violating coupling A233 approach a non-trivial fixed point, and the r- Yukawa coupling 
approaches a trivial fixed point. All other possible fixed points are either unphysical or 
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unstable in the infrared region. 

It is well known that the minimal supersymmetric standard model suffers from the so- 
called [L problem associated with the bilinear term connecting the two Higgs doublet super- 
fields in the superpotential. A simple solution to this problem is to postulate the existence 
of a Higgs singlet superfield, and to couple it to the two Higgs doublets in the superpotential 
via a dimensionless trilinear coupling. When the Higgs singlet obtains a vacuum expectation 
value, a bilinear term involving the two Higgs doublets is naturally generated [12]. However, 
this leads to additional trilinear superpotential couplings in the model, the so called non- 
minimal supersymmetric standard model (NMSSM). Furthermore, if we do not postulate 
R p conservation, then there is also an additional lepton number violating superpotential 
coupling in this model as compared to the MSSM with baryon and lepton number violation. 
Because of these additional trilinear couplings, it is important to study the infra-red fixed 
point structure of the NMSSM, and analyze the effect of these additional Yukawa couplings 
on the infrared behavior of the other Yukawa , and the baryon- and lepton-number violating 
couplings, and contrast it with the situation that obtains in the MSSM. Some preliminary 



studies of the NMSSM with B and L violation were carried out in ||13|| . In this paper we 
carry out a detailed study of the renormalization group evolution of the Yukawa couplings 
of NMSSM, including the B and L violating couplings. We shall include all the third gen- 
eration Yukawa couplings, as well as the highest generation B and L violating couplings 
in our study, and analyze the situation where all of them could simultaneously approach 
infrared fixed points. We shall investigate both the true infrared fixed points, as well as 
quasi-fixed-points of these couplings. In particular, we shall carry out a detailed stability 
analysis of the infrared fixed points of these couplings. Furthermore, corresponding to the 
B and L (and R p ) violating Yukawa couplings of the NMSSM, there are soft supersymmetry 
breaking trilinear couplings (the A parameters) whose renormalization group evolution and 
infrared fixed point structure has not been studied so far. We shall, therefore, also study 
the renormalization group evolution of these soft supersymmetry breaking A parameters, 
including those corresponding to the third generation Yukawa couplings, and obtain the 
simultaneous infrared fixed points for them. 

The plan of the paper is as follows. In Sec. II we describe the model and write down all 
the dimension four and five baryon and lepton number violating couplings in the NMSSM. 
We then derive the renormalization group equations (RGE's) of interest to us, which include 
the equations for the dimension four baryon and lepton number violating Yukawa couplings. 
We then carry out a detailed analytical study of the true infrared fixed points of the Yukawa 
couplings in full generality. Within the context of grand unified theories, one is led to the 
situation where B and L violating Yukawa couplings may be related at the GUT scale, and 
one may no longer be able to set one or the other arbitarily to zero. We, therefore, initially 
include both baryon and lepton number violating couplings in our RG equations. The fixed 
point analysis of such a system of RG equations leads to the crucial result that the only 
stable fixed point is the one with simultaneous non-trivial fixed point values for the top- 
and bottom-quark Yukawa couplings and the B- violating coupling A233, and a trivial one for 
all other couplings. Thus, non-trivial simultaneous fixed points for the B and L violating 
Yukawa couplings are ruled out by our analysis of the NMSSM. This is result is analogous to 



the corresponding result obtained in the MSSM ||11|| . We then study the fixed points of the 



corresponding soft supersymmetry breaking trilinear couplings of this model. In Sec. Ill we 
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algebraically study the simultaneous quasi-fixed points of all the third generation Yukawa 
couplings of the minimal supersymmetric standard model with B violation, as well as those 
of the corresponding soft supersymmetry breaking trilinear couplings. Since the quasi-fixed 
point limit is formally defined as the Landau pole of the Yukawa coupling at the GUT scale, 
it provides an upper bound on the corresponding Yukawa coupling. In Sec. IV we present 
the numerical results for the renormalization group evolution and the quasi-fixed points for 
the nonminimal supersymmetric standard model with B violation. In Sec. V we summarize 
our results and present the conclusions. 



II. RENORMALIZATION GROUP EQUATIONS AND INFRA-RED FIXED 

POINTS 

A. Baryon and lepton number violation in NMSSM 

In this section we study the true infra-red fixed points of the Yukawa couplings and the A 
parameters of the NMSSM with B and L violation. We begin by recalling the basic features 
of the NMSSM with baryon and lepton number violation. The superpotential of the model 
is written as 

W = (huUQllf^ + (h D ) ai Q0 R H 1 + {h E ) ah L a I E h R H 1 + XNH X H 2 - ^N 3 . (1) 

where L, Q, E, D, U denote the lepton and quark doublets, and anti-lepon singlet, d-type 
anti-quark singlet and u-type anti-quark singlet, respectively. In Eq. ([!]), (hu) a b, (^d)o6 
and (h,E)ab are the Yukawa coupling matrices, with a, b, c as the generation indices. Gauge 
invariance, supersymmetry and renormalizability allow the addition of the following L and 
B violating terms to the superpotential (H): 

W L = \ a NL a H 2 + l -\ abc LlL b L T R + X' abc L a L Q b L D R , (2) 
W B = \K bc D a n^ R U R , (3) 



where the notation |T3 is standard. We note that there is a additional L-violating term 



with the dimensionless Yukawa coupling A a in (g) which does not have an analogue in the 
MSSM. This term can be rotated away into the R-parity conserving term \NHiH 2 via 
an SU(4) rotation between the superfields Hi and L a . However, this rotation must be 
performed at some energy scale, and the term is regenerated through the renormalization 
group equations. The Yukawa couplings \ a b c and \" bc are antisymmetric in their first two 
indices due to SU(2) L and SU(3)c group symmetries, respectively. 

The dimension-4 terms in the supertentials (0) and (0) are the most dangerous terms for 



nucleon decay, and some of them must be suppressed. This leads to constraints |14| on the 
different couplings X a bc, A^ fec , and A" fec , but considerable freedom remains for the various B 
and L violating couplings. Furthermore, there are dimension-5 operators which may lead to 
nucleon decay suppressed by 1/M, where M is some large mass scale at which the B and L 
violation beyond that of NMSSM (and MSSM) comes into play. Some of these dimension-5 
operators may also lead to unacceptable nucleon decay if their coefficients are of order unity, 
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and therefore must be suppressed. We tabulate here all the dimension-5 operators which 
are allowed by the SU(3)c x SU{2) L x U(l)y gauge symmetry and the particle content of 
the NMSSM. These are 



01 = 


[QQQL] F , 


02 = 


[UUDE] F , 


03 = 


[QQQH X \ F , 


4 = 


[QUEHxIf, 


05 = 


[LLH 2 H 2 ] F , 


06 = 


\LExEL 2 E 2 \ F 


07 = 


[H 2 H 2 E*] D , 


08 = 


[H 2 H X E\ D , 


09 = 


[QUL*] D , 


010 = 


- [UD*E} D , 


011 = 


-- [LH 2 NN] F , 


012 = 


-- [LLEN] F , 


013 = 


- [LQDN] F , 


014 = 


= [DDUN) F , 


015 = 


- [LH 2 N*] D , 







(4) 



where we have suppressed the gauge and family indices. We note that the baryon and lep- 

ton number violating operators 0i, , O w are the same as the corresponding operators 

in MSSM and are subject to constraints similar to those in MSSM, whereas the opera- 
tors 0n, ,0i5 are additional dimension-5 B and L violating operators specific to the 

NMSSM. We note that in NMSSM the simutaneous presence of the combinations LQDN 
and DDU, or LQD and DDUN needs to be forbidden, since these could lead to fast proton 
decay. We do not consider the dimension-5 operators (f|) any further in this paper, and 
restrict our attention only to the dimension-4 terms in (^), and (||), respectively. 

Corresponding to the terms in the superpotentials ([U), (||) and (||), there are the soft 
supersymmetry breaking trilinear terms which can be written as 



soft 



(AuUihuUQiURHz + {A D ) ab {h D ) ab Q a L D R H x 



+ (AEUihEULlE^ + A x XNHiH 2 - ^iV 3 



+ 
+ 



~ ~ j ~ ~ ~ Q ^_ ^_ ■— (I 

{A~ x ) a KNL a L H 2 + -(Ax)abcKbcL a L L L E R + (A y ) abc X' abc L a L Q L D R 

I ~ a ~ b ~ c 

-(Ax")abcKbcD R D R U R 



(5) 



where a tilde over a matter chiral superfield denotes its scalar component, and the notation 
for the scalar component of the Higgs superfield is the same as that of the corresponding 
superfield. In addition there are soft supersymmetry breaking gaugino mass terms with the 
masses M i; with i — 1,2,3 corresponding to the gauge groups £/(l)y, SU(2)l, and SU(3)c, 
respectively. 

The third generation Yukawa couplings are the dominant couplings in the superpotential 
([!]). Therefore, it is natural to retain only the elements (hu)33 = h t , (h D ) 33 = h b , (h L ) 33 = 
h T in each of the Yukawa couplings matrices hjj, hp, h^, setting all other elements equal 
to zero. Furthermore, there are 39 independent L violating trilinear couplings A a , X abc and 
X' abc in (Q). Similarly, there are 9 independent B violating couplings X" bc in the baryon 
number violating superpotential (||). Thus, we would have to consider 44 coupled nonlinear 
evolution equations for the L violating case and 14 coupled nonlinear equations for the B 
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violating case, respectively to study the renormalization group evolution of these couplings 
in the NMSSM. It is clear that there is a need for a radical simplification of these equations 
before we can think of studying the evolution of the Yukawa couplings in the NMSSM with 
B and L violation. 

In order to render the Yukawa couping evolution equations tractable, we, therefore, need 
to make certain plausible assumptions. Motivated by the generational hierarchy of the 
conventional Yukawa couplings, we shall assume that an analogous hierarchy amongst the 
different generations of B and L violating couplings exists. Thus, we shall retain only the 
couplings A3, A 233 , A333, A233, and neglect the rest. We note that B and L violating cou- 
plings to higher generations evolve more strongly because of larger Higgs couplings in their 
evolution equations, and hence could take larger values than the corresponding couplings 
to the lighter generations. We also note that the experimental upper limits are stronger for 
the B and L violating couplings in MSSM with lower indices ||14|| . 



B. Renormalization group equations 

We are interested in the one-loop renormalization group equations for the dimension- 
less trilinear Yukawa couplings in the superpotential ([]]), (^|) and (|3]). For a general N=l 
super symmetric theory with a trilinear superpotential term / a b c $ a $ 6 $ c involving chiral su- 
perfields $ a , $ b , $ c , the evolution of the couplings f a bc with the scale parameter \i is given 



by the RGEs @ 



167T 2 |^ = Yafebc + Tf/oac + YJabe, (6) 

where 7^ are the elements of the anomalous dimension matrix, and sum over repeated indices 
is understood. The anomalous dimensions are given by 

i: = \Y.f^f ebc -25:g\Ct (7) 

b,c 

to one loop order. The sum over A represents a sum over all dominant gauge couplings, 
and is the quadratic Casimir of the representation of $ a under the gauge group with 
coupling g A : 

(T£ Ti) b a = C£ 8 b a . (8) 

Here, T£ is a matrix in the R representation for the group labelled by A. Pictorially, the 
RG evolution of the trilinear coupling can be described as an insertion of the anomalous 
dimension correction on each external leg. We have calculated the anomalous dimension for 
the various superfields for the NMSSM with baryon- and lepton-number violating couplings. 
These are summarized in Table |I[ The renormalization group equation for the Yukawa cou- 
plings hu,h,D,hE of the superpotential (JTJ) are obtained from (g) with the index c belonging 
to a Higgs field. The general form of the RGEs are 
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$ a ' b NMSSM L violation B violation 



N,N 4A 2 + 4fc 2 A*Aj 

Li, Hi - X ah (hE)ab + 3\' iab (hD)ab + AA* 

hj h L rf L -lgl-^gl hab^ jab + 3\' iab \' jab + AjA^j- 

^ 2^ E -| 5l 2 A obi A o6 ,- 

f) c 9ft t 8_2 2 2 9 \'a6i\' 9\"ia6\" 

^ D ^D - 3 #3 - 15 fl 2A A abj ZA A jab 

% l^hv-lgl-yi - A'^A;, 

Qi,i huhjj + h D h], - lg 2 3 - \g\ - ±g\ \' aib \' ajb 

HuHt Tr{h E h ] E + Zh D lJ D ) + \ 2 -lg 2 2 -^gl 

H 2 ,H 2 3Tr(huh\j) + A 2 - \g\ - ^g\ A,V 



TABLE I. The anomalous dimensions 7?^ in the non-minimal supersymmetric standard model 
with lepton and baryon number violating couplings. Here a, b are flavour indices. 

16vr 2 ^(Ma6 = (hu) lbl 2: + {hu)ail% + (huU-rgl (9) 

16n2 d]Ljl {hD)ab = (Mi67 £ + (^)-^S + ( D Ulm + 4&7fe, (10) 
167r 2 ^-(/i L ) ab = (h L ) ib j% + (h L ) a all + (HlUihI + AiaftTffa- ( n ) 

The evolution of gauge couplings ^ (i = 1,2,3 denoting the C/(l)y, SU{2) L and SU(3)c 
gauge groups, with GUT normalization for i7(l)y gauge group) in NMSSM with B- and 
L-violation is same as in MSSM, as this evolution is unaffected at the one-loop level by 
the presence of a chiral singlet superfield, or B- and L-violating couplings. These evolution 
equations are 

l 67r 2^ =& 3 ; i=123 (12 ) 

aln/z 

where &, are the beta functions for the respective gauge couplings with b\ = 33/5, b 2 = 
1, 63 = —3. The corresponding one- loop renormalization group equations for the gaugino 
masses M i; i = 1, 2, 3 can be written as 

16. 2 ^ = 2 ?l 2 W, z = l,2,3. (13) 
a In// 

Retaining only the third generation Yukawa couplings and the highest generation baryon 
and lepton number violating couplings, the RGEs for the Yukawa couplings and the R-parity 



violating couplings in the NMSSM can be written as |L6 
dh t h t 



(iln/i 167r 2 



2 + hi + A 2 + A 2 + A333 + 2A 2 2 3 



16 



-9s 



^ + ^ 2 



dhh 



1 



din fi 167r 2 



h 2 t +6h 2 b + h 2 T + \ 2 + 6\' 3 2 33 + 2\" 2 



233 



ft, 6 + AA 3 A 



dln/j 167T 2 



16 2 o 2 7 2 



3/i 2 + 4h 2 T + A 2 + A 2 



4A 2 33 + 3A' 2 



338 



9 



- ( 3^ 2 + 5^ 



dX 



1 



3fr 2 + 3h 2 b + h 2 T + 4A 2 + Ik 2 + 4A 2 ) A + 3/i 6 A 3 A 



6^3^333 



o 2 i 2 



A 



dn 
(iln/i 

c*A 3 
(iln/i 



ft 



167T 2 
1 



6A 2 + 6k 2 + 6A1 



dX 



233 



167T 2 

-|-3/ifc AA333 
A233 



3h 2 + h 2 T + 4A 2 + 2k 2 + 4A 2 + A 
3 



2 

233 



+ 3A333 



3<fe 2 



A s 



(iln/i 

^333 
(iln/i 



dA 



233 



167T 2 
1 

167T 2 

+/i b AA 3 
\" 

A 233 



4/ i 2 + A 2 + 4A 2 33 + 3A 3 2 3 - 
h 2 + 6h 2 + h 2 T +~Xl + A 2 33 



3<fe 2 



16 



7 



J93 + 3^2 + J^l 



+ 6A 



x 333 



'2 
333 



2A 



2h; + 2hi + 2A 



'2 

333 



6 A 233 



8^ 3 2 



233 



4 2 
5^ 



A 



333 



(14) 

(15) 

(16) 

(17) 
(18) 



(19) 
(20) 



(21) 



(22) 



(iln/i 167r 2 

We note that since the difference between the one-loop and two-loop results jT7| for the 
infrared fixed points in MSSM is less than 10%, we shall use one-loop renormalization group 
equations in the study of infrared fixed points in the NMSSM in this paper. 

We now come to the evolution equations for the soft supersymmetry breaking trilinear 
parameters in the potential (|5|). The one- loop RGEs for these parameters can be deduced 
from the general expressions in ref. QlS | . In this paper we shall assume the same kind of 
generational hierarchy for these trilinear parameters as was assumed for the corresponding 
Yukawa couplings. Thus, we shall consider, besides A\ and A K , only the highest generation 



trilinear coulings (A u ) 33 = A t , (A D ) 33 = A b , (A L ) 33 = A T , (A~ x ) 3 = A~ x (A 



UJ233 



= A 



A233 ' 



(Ay) 333 = Ay , (Ay/) 233 = A\n , setting all other elements equal to zero. With this 

333 233 

assumption the RGEs for the soft supersymmetry breaking trilinear parameters can be 
written as 

(1A ' 1 ( 1 ,:? " ' A,X 2 + A~~Xl + A x , 



(iln/i 8 



7T^ 



6A t hf + A b h 



^3^3 



A 

A 333 333 



+ 2Ay> X 



112 

233 



X6 1 3 \ 

-—M 3 g 2 3 - 3M 2 g\ - -M l9 2 , 

^ (^/i 2 + 6A b h 2 b + A T h 2 T + A\\ 2 + ^A b \f 33 + ^Ay 333 X 

+2A A" 2 - ^ fcA ^ A 333 , 4\AA 3 A'333 ^A^AAsA^a 

A * 33 233 2/i 6 h b 2h b 
~M s g* - 3M 2 g\ - ^M 1 ^j , 

1 ^ 3A b h 2 b + 4A r /> 2 + A A A 2 + ^A 2 + + 7Aa2mA 



/2 

'333 



8tt 2 

+3A Ak3 A 3 2 33 - 3M 2 g\ - ^M l9 f 

3A t h 2 + 3A b h\ + A r /i 2 + 4A A A 2 + + 2A K /t 2 
7 ^a 3 -^3 3Ab/i b A 3 A 3 33 _ 3A A /ibA3A 33 3 3A A3 /i fe A 3 A 33 3 



» 2 

A 2 33 A 233 



8tt 2 
-f 

-3M 2 g 2 2 - ±M l9 \ 

(a x X 2 + A k k 2 + A^lf) , 



A 2A 2A 

3., o 



6 



i / 74 \2 /I- \2 

8^ (s^ 2 + A r /. 2 + ZM- + 2A K . 2 + - + 4^ 3 A 2 

4 x2 , OA W2 | 3A A /l b AA 3 33 3^ 3 /l 6 AA 33 3 
+^A 233 ^233 + ^A^Agaa + 

3 Ay /ifoAAooo „ 3 Q \ 

+ ^333 333 _ 3M ^2 _ Mi ^2 ) 

A 3 5 / 

^2 + ^3 + ^ + ^A 233 AL + 3A A , 33 A' 2 3 



3M 2 g\ - -M x9 \ 



5 

1 

-f 



A^ 2 + + ^U 2 + A, 3 A 2 + A A233 A 2 33 

+ QAy A 3 2 3 + 2A A » A- + Abk f X3 + 

O A 333 333 -^233 233 o \/ \/ 

z ^ A 333 A 333 

A x , h b XX 3 16 2 2 7 . . 2 \ 

^, y ^ 3 " 3M 2 (? 2 2 - — M^ 2 j , 

L (2A t /i 2 + 2A 6 /i 2 + 2A A ,A 3 2 33 + 6^A'«3 
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-8M 3 g 2 3 - \m i9 1\ . (31) 

Given the evolution equations (|14D - fl22|) for the Yukawa couplings and the evolution 
equations (|23|) - ( pT|) for the A parameters, we are now ready to study the RG evolution and 
infra-red fixed points of the NMSSM with B and L violation. 



C. Infrared fixed points for Yukawa couplings 

In this section we consider the infrared fixed points for the Yukawa couplings and the 



baryon and lepton number violating couplings of the NMSSM. The infra-red fixed points JT3 
of the renormalization group equations ([14] ) - (22) have been studied in the limit of ignoring 
the r- Yukawa coupling h T , and by considering either the baryon number violating Yukawa 
coupling A233, or the lepton number violating Yukawa couplings A3, A233 and A 333 . In the 
analysis that follows, we shall consider the evolution equations for h t and hb together with 
the evolution equation for h T . Furthermore, we shall also entertain the possibility of si- 
multaneous presence of B and L violating couplings in the renormalization group equations 



(|14]) - (122|). We do this in order to investigate as to whether such a system of equations does 
have acceptable infrared fixed points. However, without loss of generality, we shall assume 
that there is a hierarchy of lepton number violating couplings, and consider only one lepton 
number violating coupling, together with the baryon number violating coupling A 233 , at a 
time. Thus, we shall consider three different cases, i.e., we shall take A 3 ^> A233,A 333 , or 
A233 3> A 3 , A 333 , or A 333 3> A 3 , A233. In order to study the infrared fixed points, we define the 



:«3 



following ratios of the dimensionless trilinear couplings and the SU(3)c gauge coupling: 

(32) 

(33) 

We shall first consider the RG evolution of Yukawa couplings in the superpotential (|l]), 
with B violation arising from (|3|), and L violation arising from the first term in (|2|). With 
the definitions (|32|) and (0), and retaining only the SU(3)c gauge coupling, the one-loop 
renormalization group equations for h t , 
form 





hi 




h 2 




k * R 
—5-, -Ka — 


A 2 R 


K 2 


Rt = 


Rb = 


2 ' 


R T = 








9V 




9i 




9i 


9i 




R 3 = 


\ 2 

A 3 


R = 


X 2 

A 233 
2 ' 

9i 


R' = 


-^333 d" 
2 ' K 
91 


\" 2 
A 233 





hb, h T , A, k , and A3, A 233 can be written in the 



dR+ 



d{ 


— In /i 2 ) 




dRb 


d{ 


— In /1 2 ) 




dR T 


d{ 


— In /1 2 ) 




dR x 


d{ 


— In /i 2 ) 




dR K 



16 



— + h ) -GRt-Rb-Rx-Rs- 2R 



d(— In jj? 



a 3 R t 
a 3 R b 

a 3 R T [b 3 -3R b -4R T -R x -R 3 }, 



16 

y + & 3 ) -Ri 0/?,., - R T - R x -2R 



a 3 R x [b 3 - 3R t - 3R b - R T - 4i? A - 2R K - 4i? 3 ], 
a 3 R K [b 3 - 6R X - 6R K - 6^3], 



(34) 
(35) 
(36) 
(37) 
(38) 
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dR 3 



d( — In /i 2 ) 

dR" 
d(-ln/z 2 ) 



a 3 R 3 [b 3 - 3R t R T 4i? A - 2i? K - 4i? 3 ] 
a 3 R" [(8 + 6 3 ) - 2Rt - 2R b - 6if], 



(39) 
(40) 



where 63 = -3 is the beta function for g 3 in NMSSM (or MSSM), and a 3 = ^/(16tt 2 ). 
Choosing the basis of the ratios as Ri = (R tl Rb, R T , R\, R K , R , R 3 ), we can rewrite the RG 
equations (|3~4D - (|40|) in the form {t = — In /1 2 ) 



dRi 
~dt 



a 3 Ri[{ri + b 3 ) - ^ S'y-Rj], 



(41) 



where r» = J2k 2Cfl, Ci? is the QCD quadratic Casimir for the various fields {Cq = C\jc = 
Cd c = 4/3) and the sum is over the representation of the three fields associated with the 
trilinear coupling that enters the definition of Ri, and S is a matrix whose entries are the 
numerical coefficients (the wave function anomalous dimensions) of R^s in the evolution 
equations (|3~4]) - (|40|). A fixed point is, then, reached when the right hand side of Eq. ( f41~D 
is for all i. If we were to write the fixed point solutions as R*, then there are two fixed 
point values for each coupling: R* — 0, or 



{ri + b^-Y.Si-jR* =0. 

j 

It follows that the non-trivial fixed point solution is 

r: = Y,( s ~ 1 u^+^)- 



(42) 



(43) 



The anomalous dimension matrix S that enters Eq. (|4*3"|), which we denote by Sbli in this 
case, is easily seen to be 



Inverting the matrix 
solution: 

RI 
R"* 



S 



BLl 



6 10 10 2 1 

16 110 2 

3 4 1 1 

3 3 1 4 2 4 

6 6 6 

2 2 6 

3 1 4 2 4 



(44) 



and substituting in Eq.(EJ), we get the following fixed point 



29 

76' 
161 

228' 



T7* 



Rt 



0, R* 

44 
~19' 



31 
76' 



RI 



18 
19' 



7?* 



33 
38' 



(45) 
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We note that R*, R% < 0, and, therefore, this fixed point solution is physically unacceptable. 
We conclude that a simultaneous fixed point for the B and L violating couplings A 2 3 3 and 
A 3 , and the Yukawa couplings h T , h b , h t ,X,K, does not exist. 

Next we consider the RG evolution of Yukawa couplings in the superpotential ([!]), with B 
violation arising from ([5[), and L violation arising from the second term in @. The one- loop 
renormalization group equations for h t , hb, h T , A, k , and A233 and A 233 can be written in 



the form ([41]) , with the anomalous dimension matrix S given by 



S 



BL2 



6 10 10 2 

16 110 2 

3 4 1 4 

3 3 1 4 2 

6 6 

4 4 

2 2 6 



(46) 



where the ordering of the ratios is Ri = (R t , Rb, R T , R\, R K , R, R")- From 
get the fixed point values 



and we 



R* t 



R* 



R 



3 

6 "4 



5 9 7 

p* p* ' p* 

K T - p X x - -, K K - -, 



5 

12' 



(47) 



Since R*, R* x < 0, this fixed point must also be rejected as being unphysical. Thus, a 
simultaneous fixed point for the B and L violating couplings A233 and A233, and the Yukawa 
couplings h T , h b , h t , A, k, does not exist. 

Finally, we consider the RG evolution of Yukawa couplings in the superpotential ([[]), 
with B violation arising from (|3|), and L violation arising from the third term in (|2]). The 
one-loop renormalization group equations for h t , hb, h T , A, n , and A 333 and A 233 can again 
be written in the form (^TJ), with the anomalous dimension matrix S given by 



S 



BL3 



6 10 10 12 

16 110 6 2 

3 4 1 3 

3 3 1 4 2 

6 6 

1 6 1 6 2 

2 2 2 6 



(48) 



where the ordering of the ratios is Ri = (R t , Rb, R T , R\, R K , R R ")- This results in the fixed 
point values 



r' * 



292' b 

™, R"< 
219' 



98 



219' 
611 

876' 



RI 



285 
292 



> R\ = 0) R* K 



1 

2' 



(49) 
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which must also be rejected as being unphysical. We conclude that B and L violating 
couplings of the highest generation cannot simultaneously approach a non-trivial fixed point 
in the NMSSM. This is one of the most important conclusions that we draw from the analysis 
of the renormalization group equations of NMSSM, and is analogous to the corresponding 
result in MSSM 0. 

Given this result, it is now natural to consider the possibility of having either B or 
L violation, but not both simultaneously, involving the trilinear couplings with highest 
generation indices in the RG evolution in the NMSSM. 



1. Fixed points with baryon number violation 

In this section we shall consider the infrared fixed points of the NMSSM with B viola- 
tion. Thus, we shall consider the Yukawa couplings h t , hb, h T , A, k, and the baryon number 
violating coupling A 233 . We order the ratios of the trilinear couplings to the gauge coupling 
g 3 as Ri = (Rt, Rb, R T , R\, Rk, R ), resulting in the anomalous dimension matrix 



S 



Bl 



6 10 10 2 

16 110 2 

3 4 1 

3 3 1 4 2 

6 6 

2 2 6 



(50) 



leading to the fixed point values 



Ri 



9 

16' 



7?* 



11 

16' 



Rl 



Rl 



41 
16' 



T?* 



33 
16' 



R"* 



5 

12' 



(51) 



which must be rejected as being unphysical. We are, therefore, constrained to consider fixed 
point with one of the couplings approaching a zero fixed point value, with all others attaining 
a non-trivial fixed point value. We first consider a fixed point with R* = 0, with all others 
obtaining a non-zero fixed point value. In this case we obtain the fixed point values 



Rl 



0, R$ 



Rl = - 



23 



r: 



19 



R 



5 

12' 



(52) 



which is a physically unacceptable fixed point solution. Proceeding in this manner, it can 
be shown that there is no acceptable infrared fixed point solution with one of the couplings 
ht-, hb, h T , A, n, A 233 approaching a zero fixed point value, and the remaining approaching a 
non-zero fixed point value. 

We are, therefore, led to the consideration of a fixed point in which two of the Yukawa 
couplings approach a zero fixed point value, with all others attaining a non-zero fixed point 
value. We try the fixed point with R* x = R* K = 0, with the rest approaching a non-zero fixed 
point value. In this case we obtain the fixed point values 



31 
292' 



^~ 73' 



Rl 



285 
292' 



R 



611 

876' 



(53) 
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which must be rejected as a fixed point solution. Taking other pairs of Yukawa couplings to 
have a zero fixed point values, with the remaining obtaining non-zero values, we find that 
there is no acceptable fixed point solution in this case. 

We next consider the situation in which three of the couplings attain zero fixed point 
values, with the remaining attaining a non-zero fixed point value. Taking R* T = R\ = R* K = 0, 
we obtain the fixed point 



p* r>* r> * 



17' 



17' 



77 
102' 



(54) 



which is a physically acceptable fixed point solution. We find that there is no other accept- 
able fixed point solution with three of the couplings obtaining fixed point values, with the 
remaining approaching a non-trivial fixed point value. 

We now consider the case when four of the trilinear Yukawa couplings approach a trivial 
fixed point, with the remaining two having a non-trivial fixed point value. In this case we 
find the following acceptable fixed point solutions: 



p* p* p* p 

It^ — lt\ — it., — it 



0, B* 



1 



1 

3' 



and 



r: 



r: 



p* 



o, r; = ~ r"* = H 

' * 8' 24' 



(55) 



(56) 



Since there are more than one theoretically acceptable IRFPs in this case, it is necessary 
to determine, which, if any, of these fixed points is more likely to be realized in nature. To 
this end, we must examine the stability of each of the fixed point solutions fl54]), (B3) and 



The infra-red stability of a fixed point is determined by the sign of the quantitites 



* = r 



SijR^-in + h) 

j =771+1 



(57) 



for those couplings which have a fixed point value zero, R* = 0, i = 1,2, ...,m, and by the 
sign of the eigenvalues of the matrix 



.4 



—RiSij, i 



m + 1, 



(58) 



for those couplings which have a non-trivial IRFP ||, where R* is the set of the non- 
trivial fixed point values of the Yukawa couplings under consideration, and is the matrix 



appearing in the corresponding renormalization group equations fl4"i ) for the ratios R{. For 
stability, we require all the Aj's to have negative sign, and the eigenvalues of the matrix (|58"D 
to have negative real parts (note that the QCD /3-function 63 is negative). For the infra-red 
fixed point ([54]) , we find from Eq. 



(59) 



A - 19 A - 21 A 
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corresponding to R* T = 0, R\ = 0, and R* = 0, respectively and for the eigenvalues of the 
matrix ( [58f) 

x 10 no ^ -273 - ^43113 -273 + Vim 

= _ 51 ~ -°' 2 ' = 306 " - 1 " 6 ' Ae = 306 * -°' 2 ' (60) 

corresponding to the non-trivial fixed point values for R"*, R£ and R£, respectively. Since 
all the Aj, i = 1,2, 3, 4, 5, 6 are negative, the fixed point ([54]) is infra-red stable. 

Next we consider the stability of the fixed point flSHp. Since in this case R\ = R* K = 
R* = R"* = 0, we have to obtain the behaviour of these couplings around the origin. This 
behaviour is determined by the quantities ( |57| ) which, in this case, are 

Ai = ~ A 2 = -~ A 3 = -1 A 4 = H (61) 

thereby indicating that the fixed point is unstable in the infra-red region. For completeness 
we also obtain the behaviour of Rb and Rt around their respective fixed points governed 
by the corresponding eigenvalues of the matrix ([58]). We obtain for the eigenvalues of this 
matrix 

A 5 = A 6 = --. (62) 

Although A5, \§ are negative, the fact that A4 > implies that the fixed point (|55|) is 
unstable in the infra-red region. Thus, this infra-red fixed point with trivial fixed point 
value for the baryon number violating coupling A'233 will never be realized at low energies, 
and must be rejected. Similarly, it is straightforward to see that the fixed point ( |56"D is 
unstable in the infra-red, and must be rejected. 

One may also consider the case where the couplings A233, k, A, h T , attain trivial fixed 
point values, whereas h t attains a non-trivial fixed point value. In this case we find Rl = 
7/18, the well-known Pendleton-Ross type fixed point p. The stability of this fixed point 
solution is obtained by simply considering the quantities fl57|) 



K=rK S BiU R Z-( r i + h)}, ^ = 1,2,3,4,5, (63) 

where the matrix Sbi is the one which appears in (pD[ ) but with the ordering of the ratios 
as Ri = (R T , R\, R K , Rb, R , Rt), which yields 

25 , 35 , 38 

A! = -1, A 2 = --, A 3 = -1, A 3 = -, A 5 = -, (64) 

thereby rendering this fixed point unstable. 

Finally, one may consider the case where R"* = 0, with Rt, R b , R T , R\, R K , attaining 
non-trivial fixed point values. This is the case of the NMSSM with R p conservation and 
with all the third generation Yukawa couplings taken into account. In this case, we find the 
fixed point solution: 

d* _ 37 33 5 51 43 

Rt - ^ , R b - T% R T - - g, R\ - - jg, R k - ^, (65) 
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which must be rejected as being unphysical. 

Thus, we have shown that the only fixed point which is stable in the infra red region is the 
baryon number, and R p , violating solution (|54]). We note that the value of R$ corresponding 
to this solution is lower than the fixed point value of 7/18 for NMSSM (and MSSM) with 
baryon number, and R p , conservation. 



2. Fixed points with lepton number violation 



In order to investigate the possibility of having a stable fixed point with lepton number 
violation, we now study the renormalization group evolution for the lepton number, and R p , 
violating couplings in the superpotential @. We shall consider the trilinear couplings A 3 , A 233 
and A 333 , together with the Yukawa couplings h t , h b , h T , A, k. As indicated in Sec. II, we 
shall assume, without loss of generality, that there is a hierarchy of lepton number violating 
couplings, and consider only one lepton number violating coupling at a time. Thus, we shall 
consider three different cases, i.e., we shall take either A 3 ^> A 233 , A 333 , or A 233 A 3 , A' 333 , or 
-^333 > ^3, A233, respectively. 

We first consider the case when A 3 is the dominant lepton number violating coupling. 
Ordering the ratios of the squares of the trilinear couplings to the square of the gauge 



coupling g 3 as Ri = (R t 
case is given by 



Rb, Rt, R\, Rk 



R 3 ), the anomalous dimension matrix, Slu, in this 



S 



Lll 



leading to the fixed point values 



32 

K t - — , K b - U, K T 



6 10 10 1 

16 110 

3 4 1 1 

3 3 1 4 2 4 

6 6 6 

3 1 4 2 4 



(66) 



15' A 



17 



7?* 

-n,„ 



43 
30' 



172 



45 



(67) 



Since i?*, i? 3 < 0, this fixed point is unacceptable. Thus, there is no nontrivial infrared fixed 
point for both the lepton number violating coupling A 3 , and all other trilinear couplings in 
the NMSSM. 

Next we consider we consider a trivial fixed point for one of the couplings, and nontrivial 
one for the. Taking Rl = 0, which is relevant for low values of tan/3, the anomalous 
dimension matrix is 



S 



L12 



6 10 1 

4 10 1 

3 14 2 4 

6 6 6 

3 14 2 4 



(68) 
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which is singular. This solution corresponds to a fixed line or a surface []. Continuing in this 
manner, and taking each one of the couplings to have a trivial fixed point value with the 
rest a nontrivial fixed point value, we find there is no acceptable infrared fixed point in this 
case. 

We now try to obtain a fixed point with two of the couplings approaching a trivial fixed 
point, while the rest having a non-trivial fixed point. Taking R* x — R* — 0, we obtain the 
fixed point values for the remaining couplings 



11 11 

IX, — — -, Xti — — -, It T 



24' 



24' 



7 ~ 

P* 

, -Tin 



(69) 



which is an unacceptable fixed point. Continuing in this manner we find that there are no 
acceptable infrared fixed points in this case. We have also checked that for the case with three 
or four of the couplings approaching a zero fixed point value, with the remaining having non- 
trivial fixed point values, there are no physically acceptable fixed point solutions. Finally, 
for the case with all the couplings, except the top-quark Yukawa coupling, approaching a 
trivial fixed point value, we have 



lih — li-r 



T3* D* 

It„ JXn 



o,r; = — , 

' 1 18 



(70) 



which is an acceptable fixed point. However, it is easily seen that this fixed point is unstable 
in the infrared region. We, therefore, conclude that for the case with lepton number violation 
through the coupling A3, there are no acceptable infra-red fixed points. 

If on the other hand the coupling A233 is the dominant of the lepton number couplings, 
then the anomalous dimension matrix is given by 



S 



L21 



6 10 10 

16 110 

3 4 1 4 

3 3 1 4 2 

6 6 

4 4 



(71) 



in the basis _R; = (i? t , Rb, R T , R\, R K , R), leading to the fixed point values 



1 7 23 

p* p* p* 

*t - 24 , - 24 , K 



R\ 



23 



p* 



19 



R* 



2 



(72) 



which must be rejected as a fixed point. We can try with one of the couplings, namely R% 
approaching a trivial fixed point value with the rest approaching a nontrivial fixed point. 
This case is relevant for low values of tan/?. We get 



7 19 

p* n p* p* p 

ri h — U, it, — — , ix^ — , Jx\ 

b 18 6 



0, K 



1 D* 29 

-- R* = — , 

2 12' 



(73) 



1 Here we correct an error in (jDJ 
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which is also unacceptable. We may also try a fixed point with R* = 0, and a non-trivial 
fixed point for other couplings. In this case we get 



5 5 7 

D* n D* D* D P* 



3 ir 



3 
4' 



(74) 



making this also an unacceptable fixed point. Continuing in this fashion, we find that there 
are no fixed points with one of the couplings approaching a zero fixed point value, with the 
rest having non-trivial fixed point values. 

We now try to obtain fixed points with two of the couplings having a fixed point value 
zero, and the rest having non-zero fixed point values. Taking R% = R* = 0, we obtain 



z,*_20 19 
* ~~ 27' — 9~' K 



29 r* - 3 

18 _ 4' 



(75) 



which is unphysical, and hence unacceptable. Similarly, we have checked there is no physi- 
cally acceptable fixed point with any of the two couplings approaching a trivial fixed point 
value, with the rest attaining a non-trivial fixed point value. Continuing in this fashion, we 
find that there are no acceptable fixed points in the infrared region for the lepton number 
violating coupling A233. In particular, the fixed point with five of the couplings approaching 
a trivial fixed point value 



p* p* p p* p* n 

it b — it T — rC\ — ri K — it — U, 



(76) 



and R t approaching a non-trivial fixed point value, Rl = 7/18, is unstable in the infrared 
region. 

Finally, we consider the case when the coupling A333 is the dominant of the lepton number 
violating couplings. With the ordering of the couplings Ri = (Rt, Rf,, R T , R\, R K , R'), the 
anomalous dimension matrix in this case is 



Sl3i — 



6 10 10 1 

16 110 6 

3 4 1 3 

3 3 1 4 2 

6 6 

16 10 



(77) 



leading to the fixed point values 



ttx , Iti, 

* 183 b 



107 70 _ _ 1 _ 68 

— ,R t ---,R x -0,R k ---,R __ 



(78) 



which are unacceptable. As in the other cases, we try a fixed point with one of the couplings 
approaching a zero fixed point value, and the rest non-trivial values. Taking, Rl = 0, we 
have 



R* - 121 R* — 47 R - 107 R* 
210 ' 70 ' 70~ ' 



36 
35 



,R> 



17 

42' 



(79) 



which is an unacceptable fixed point. Continuing in this manner, we find that in this case 
there are only two acceptable infrared fixed points: 



18 



and 



Rl = R* r = R x = R* K = 0, R* = I R'* = 1 (80) 



R* — R x — R* K — R' — 0, R* t — -, Rl — -. (81) 



In order to determine which of these fixed points, if any, is actually realised, we need to 
check the stability of these fixed points. We first consider the fixed point (|8~0"D . In this case 
the quantities ([57]) are calculated to be 

Ai=0, A 2 = ~ A3 = ~, A 4 = -l (82) 

corresponding to the the zero fixed point values in (j80|). From ( p2|) we conclude that the 
fixed point (^) will never be reached in the infrared region. The fixed point is either a 
saddle point or an ultra-violet istable fixed point. 



The stability of the fixed point (|8~l"D is determined in a manner analogous to that of the 
fixed point (^). We find that the fixed point (|8l|) is either a saddle point, or an ultra-violet 
stable fixed point. That the stability properties of the fixed points (|80|) and (|8T|) are identical 
is a consequence of the symmetry of the renormalization group equations (|i~4l) - (|22|) . We 
conclude that there are no non-trivial stable fixed points in the infra-red region for the L 
violating coupling A333. 

To sum up, we have found that there are no IRSFPs in the NMSSM with the highest 
generation lepton number violation. This result, together with the result on the fixed point 



with baryon number violation, shows that only the simultaneous non-trivial fixed point ( |54|) 
for the baryon number violating coupling A233, and the top- and bottom-quark Yukawa 
couplings, h t and hb, is stable in the infra-red region. This result is analogous to the result 
obtained in MSSM, and shows that the additional couplings in NMSSM have no effect on 
the infrared fixed point behavior of the Yukawa and baryon number violating couplings. 

It is appropriate to examine the implications of the value of the top-quark mass predicted 
by our fixed point analysis. From (0), it is readily seen that the fixed point value for the 
top-quark Yukawa coupling translates into a top-quark (pole) mass of about m t ~ 70 sin f3 
GeV, which is incompatible with the measured value of [|l9j of top mass, m t ~ 174 GeV, for 
any value of tan/?. It follows that the true fixed point obtained here provides only a lower 
bound on the baryon number violating coupling A233 ~ 0.97. 



D. Infrared fixed points for the trilinear soft supersymmetry breaking parameters 

Having obtained the infrared fixed point structure of the Yukawa couplings of NMSSM 
with baryon and lepton number violation, we now consider the renormalization group evolu- 
tion and the fixed point structure for the soft supersymmetry breaking trilinear parameters 
Ai. Since there is only one IRSFP in this case, we shall consider the IRFPs for the A 
parameters corresponding to this case only, i.e. for A t , A b , A T , A\, A K , and Ay^. 

Retaining only these parameters, and with the definitions Ai = Ai/M% (Ai = 
A t , Af,, A T , A\, A K , A\» 33 ), we obtain from Eq. (p3|) - Eq. (|3lD the renormalization group 

equations for the relevant Ai (neglecting the SU(2)l and U(l)y gauge couplings): 
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dA* 



d(— In n 2 ) 

dA b 
d{ — In fj, 2 ) 

dA T 
d(— In n 2 ) 

dA x 
rf(-ln/x 2 ) 

dA K 
d{- In fi 2 ) 



dA 



A,,. 



«3 



a 3 



«3 



a 3 



«3 



«3 



16 

- - (QR t - b 3 )A t - R b A b - R"A Ks , 



16 



- RtA t - (6R b - b 3 )A b - R T A T - R X A X - 2R" Ay, 



d( — In fi 2 

which can be written in a compact form 

dAi 



-3R b A b - (4R T - b 3 )A\ - R X A X 
-3R t A t - 3R b A b - R T A T - (4R A - h,)A.\ - 2/?„.,l„. 
-6i? A i A - (6R K - b 3 )A K 
8 - 2R t A t - 2R b A b - (6R" - b 3 )A x 



d(— In jj? 



«3 



(83) 
(84) 
(85) 
(86) 
(87) 
(88) 



(89) 



in the basis Ai = (A T , A\, A K , A\^ , A b , At), where = (0,0,0,8,16/3,16/3), and where K 
is a matrix whose entries are fully specified by the wave function anomalous dimensions and 
Ri. A fixed point for A4 is, then, reached when the right hand side of Eq. (|89"D vanishes for 
all i. Denoting this fixed point solution by A*, we have 



(90) 



where K* is the matrix K evaluated when Ri take their fixed point values R*. From Eq. (|83] 
- Eq. ( jS8|) we see that the matrix K at the fixed point can be written as 



K* 



4R* T - h 


Rl 








3R* b 





R* 


(4i?l - 6s) 


2R* K 





3R* b 


3R* t 





6R* X 


(QK - h) 




















(6/r* - 63) 


2R* b 


2R* t 


R* 


R* x 





2R"* 


(QRt - h) 


Rl 











2R"* 




(6R* t - 



(91) 



with R* = R* x = R* K = 0, and R"*, R b , R* t given by their fixed point values in (HI). The fixed 
points A* are given by the solution of 



A* = J2( K *~ V 



ij r ji 



(92) 



with the result 



K = ~, a* x = ~, A* K = 0, a* x ,, 33 = A; = A; = 1. 



(93) 
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We have carried out the stability analysis of the fixed point fl93|) in a manner analogous to 
that of the Yukawa couplings in the previous subsection, and find it to be infra-red stable. 



We note that the fixed point values of A T ,A\" ,At,,A t in 



are the same as in MSSM 



with baryon number violation [] I 



III. QUASI-FIXED POINTS 

The infrared fixed points that we have discussed in the previous section are the true 
IRFPs of the renormalization group equations. However, these fixed points are not reached 
in practice, the range between the large(GUT) scale and the weak scale being too small for 
the ratios to closely approach the fixed point values. In that case, the various couplings may 
be determined by the quasi-fixed point behaviour [|J, where the value of various couplings 
at the weak scale is independent of its value at the large(GUT) scale, provided the couplings 
at the unification scale are large. In this section, we shall discuss the quasi-fixed point 
behaviour of the Yukawa couplings and the A parameters of the NMSSM with B violation, 
corresponding to the true fixed points that we have obtained in the previous section. 

In order to discuss the quasi-fixed points, we shall first write down the analytical solution 
for the trilinear couplings of the NMSSM in a closed form. Since the simultaneous fixed point 
for the third generation Yukawa couplings, A, k, and the baryon number violating coupling 
A233 is stable, we shall consider the quasi-fixed points for these couplings only. For this 
purpose we define 

h 2 h 2 h 2 

* ~ 16vr 2 ' b ~ 16vr 2 ' T ~ 16vr 2 ' 1 } 

X 2 ~ K 2 

Yx = 16V' Y « = 16V' (95) 



Y" 



\"2 
A 233 

167T 2 



(96) 



Then the solution of the RG equations fll~4[) - (|i~8l) , and fl2~2"[) for the Yukawa, the trilinear 
couplings A and k, and the B violating couplings can be written in a closed form [2H] 



Yi(t) 



1 + 0*^(0) f*F(t>)dt r 



(97) 



where Mq is some large initial scale, and where Yj stands for the functions Yt , Yb , Y T , Y\, Y K , 
and Y", respectively. Analogous notation holds for the functions F,. The quantities an are 
the diagonal elements of the wave function anamolous dimension matrix, and are given by 



a ti = {6, 6, 4, 4, 6, 6}, 
and the functions Fi are given by the set of integral equations 



(98) 
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FJt) = Et ^ (99) 

(1 + 6Y 6 (0) /„' F b (t')dt>y/S(l + 4Y X (0) /„' F A (t')^') 1/4 (l + 6^(0) F''P)dt/)W 

F b (t) = s - t s = j 

(1 + 6y t (0) Jo* )df )V6(1 + 4y T (0) Jj F T (^)rff )V4(1 + 4^(0) jj F x (t>)dt>y/* 

X = , (100) 

(l + ey'^o)/*^'^)^) 1 ^' 

F T (t) = = , ^ = , , (101) 

(1 + 6Y b (0) /„' F 6 (f)df )V2(! + 4 y A (0) /* F x (t')dt')l/A 

Exit) 

(i + 6? t (o) F t (t')dt<y/ 2 (i + ey 6 (o) jj F b (t>)dt>y/ 2 (i + 4y(o) jj F r (^') 1/4 

x = 7 , (102) 

(l + 6Y K (0)f*F K (t>)dt>y/z 

p u\ E K (t) (103) 

K{> (l + AY x (0)^F x (t')dt'f/^ 



FX® 



F"(t) = = ; = , , (104) 

(l + 6Y t (0)J*F t (t')dt')^(l + 6Y b (0)J*F t (t')dt')^ 



E"(t) 

)1/3(1 

where the functions E^t) (= E t (t), E b (t), E T (t), E x (t), E K (t) and E"{t)) are given by 

3 

Ej(t) = II ( X + b m a m {Q)t) Cjm/bm , (105) 



m=l 



with 



«40) = ^T, m = l, 2, 3, (106) 

LOU 



,'13 16\ / 7 16\ , ., , 

Q m = ( — , 3, y 1 , C 6m = I — , 3, y J , C rm = I -, 3, J , (107) 



c Ai =(|, 3, Oj, c«i= (0,0,0), c A » 33 , = Q, 0, 8 ) . (108) 



The solutions for the RG equations for the gauge couplings (fT2[) and the gaugino masses (|T3[) 
are well known and will not be repeated here. We note that ( ^7p gives the exact solution 
for the Yukawa couplings, while F»'s in ( p9[) - ( |104j) should in principle be solved iteratively. 

In the regime where the Yukawa couplings Y t (0), Y b (0), Y r (0), Y A (0), Y K (0), Y"(Q) -> 
oo with their ratios fixed, it is legitimate to drop 1 in the denominators of the equations 



and (|99|) - (|104|) , so that the exact solutions for the trilinear couplings approach the 



infrared quasi-fixed-point (IRQFP) defined by 

y l QFP (t)= ,.; ',>/■;.'.',.;• coo) 
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with 

F QFP {t) = Et ® fllO) 

F QFP {t) = Eb ^ fill) 

b (f*F t QFP (t>)dt'y/^ft F Q FP ( t >)dt')V^ } 

F?FP{t) = (^^(to^^a!:^^^)!^' (112) 

pQ FP (t) = Ex ^ (113) 

A (fiF t QFP (t')dt')V 2 (fiF b (t')Q FP dV)V^ 

F * FP(t) = (ft F^\t')dt'f^ (H4) 

F " QFP{t) = (£F^(tO^)V3(£F^(tO^)V3- (U5) 

We stress here that both the dependence on the initial conditions for each Yukawa coupling 
as well as the dependence on the ratios of initial values of Yukawa couplings have completely 
dropped out of the runnings in Eqs. (|109|) and ( |110| ) - (|115|) . In other words, the quasi-fixed- 



points ( |109| ) are independent of whether the B and L violating couplings and the third 
generation Yukawa couplings are unified or not. The fact that the ratios of the various 
Yukawa couplings do not enter Eqs. ( |1U9| ) - ( |115[ ) implies that these results are valid for 



any tan (3 regime. 

One can also obtain the complete analytic solutions of the RG equations for the trilinear 
supersymmetry breaking parameters Ai in an analogous manner. The expressions for these 
are lengthy, and will not be written here. Since the solutions for the quasi-fixed points must 
be obtained either iteratively from ( 109 ), or numerically from the RG equations, we shall 



instead study the numerical solutions for the quasi-fixed points for these A parameters, as 
well as the Yukawa couplings, in the following. 



A. Quasi-fixed points for the Yukawa couplings 

In order to determine the quasi-fixed points, one can solve the RG equations for the 
various couplings numerically. We shall do so in the next section. However, before doing so 
it is instructive and useful to obtain an analytical estimate of the quasi-fixed point values 
for these couplings. As stated above, we shall study the quasi-fixed point for the couplings 
h t , hb, h T , A, k, A233, since only the fixed point fl54|) corresponding to these couplings is stable. 

Writing the RG equations for (|HD - (|T8|), and ( f2~T| ) for these couplings in terms of Y's 
defined in (0) - (pop, the existence of quasi- fixed point requires |^T] 



dY t dY b dY T dY x dY K dY" , . 

0. (116) 



dt dt dt dt dt dt 
The RG equations for Y\ and Y K lead to 
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Y X = Y K = 0, (117) 

thereby implying that the quasi-fixed point values for A and k are zero. Substituting these 
values in the remaining RG equations, we get 

Y T * = [121a! + 100(« 2 - a 3 )] ~ -0.0016, (118) 

where we have used a.\ ~ 0.017, a 2 — 0.033, a 3 ~ 0.1 at the effective supersymmetry scale, 
which we take to be 1 TeV. This solution is clearly unphysical. We are, therefore, led to set 
Y* = 0, yielding the solution 

~ 47«i + 225a 2 + 200a 3 , , 

Y* = — — — — ~ 0.0052, 119 

* 425 

~ 13«i + 225a 2 + 200a 3 , . 

Y b * = — — — — ~ 0.0052, 120 

425 

Y"* = 2 45^ + 1305,) ^ (121) 

255 ' K J 

which leads to the quasi-fixed point values for these couplings 

h* t ~ 0.91, (122) 

hi ~ 0.90, (123) 

A233 - 1-08. (124) 

We note that these quasi-fixed point values for the respective couplings are similar to the 



values obtained for them in MSSM [11 



B. Quasi-fixed points for trilinear soft supersymmetry breaking parameters 

We now turn our attention to the renormalization group equations ( |23"D - ( |3~T| ) for the A 
parameters, and their quasi-fixed points. Since the quasi-fixed point values for h T , A, k are 
trivial, we cannot determine the quasi-fixed point values for A T , A\, and A K , and, therefore, 
ignore the coresponding renormalization group equations (p5|) - (0). In rest of the A 
equations, we substitute h T = A = k = 0, and obtain the equations that the remaining A 
parameters must satisfy in order for us to determine the quasi-fixed point solution: 

X 6 13 

6Y t A t + Y b A b + 2Y"Ax" - —a 3 M 3 - 3a 2 M 2 - — 5i = 0, (125) 

3 15 

16 7 

Y t A t + 6Y b A b + 2Y"A X u - —a 3 M 3 - 3a 2 M 2 - —ot x = 0, (126) 

2Y t A t + 2Y b A b + 6Y"A X „ - 8& 3 M 3 - = 0. (127) 
These yield the following quasi-fixed point solution: 
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= 47^+2250,^ + 200^, = (12g) 



4251? 

13a 1 M 1 + 225a 2 M 2 + 200£ 3 M 3 _ fi 
Al = ~ ~ 0.78ms, (129) 

2 (lliaiMi - 45a 2 M 2 + 130a 3 M 3 ) , , 

A* XII = — — ~ 1.02m,, 130 

255 Y"* 9 

where rrig is the gluino mass (= M 3 ) defined at the weak scale. For the numerical estimates 
in ( |128j ) - ( |130| ), we have used the fact that the gaugino masses scale as the square of the 
gauge couplings, and that «g — 0.041 at the grand unified scale Mq ~ 10 16 GeV. One 
must compare these quasi- fixed point values with the true fixed-point values (|93| ) . Since the 
quasi-fixed points for the Yukawa couplings represent their upper bounds, the quasi-fixed 
point values ( |128| ) - ( |130| ) provide a lower bound on the corresponding A parameters. We 
note that inspite of the presence of additional trilinear couplings in the superpotential of the 
NMSSM, the quasi-fixed point behaviour is similar to the corresponding quasi-fixed point 
behaviour in MSSM with baryon and lepton number violation. 



IV. NUMERICAL RESULTS AND DISCUSSION 



In the previous section we have obtained the approximate quasi-fixed point values for the 
Yukawa couplings and the A parameters by an algebraic solution of the corresponding RG 
equations. The RG equations are a set of coupled first-order differential equations that must 
be solved numerically to obtain accurate values for the fixed points. We have numerically 
solved the RG equations for the Yukawa couplings, and the A parameters. We now present 
the results of such a numerical analysis. 

In Fig.l we show the fixed point behaviour of the top-quark Yukawa coupling as a function 
of the logarithm of the scale parameter fi. We have included the evolution equations for 
the b-quark, the r-lepton Yukawa coupling, and the trilinear couplings A and k, as well as 
the B-violating coupling A 233 , and evolved them from the initial value ln(n/Mo) = to the 
final value Iti(/j,/Mg) — —33 at the weak scale, for the numerical solution. It is seen that 
for all h t ~ 1 at the GUT scale, the top-quark Yukawa coupling approaches its quasi-fixed 
point at the weak scale. We note that the numerical evolution of fixed point approaches 
but does not exactly reproduce the approximate analytical value in ( |122|) . In Figs. 2 and 3 
we present the corresponding approach to the quasi-infrared fixed point for the couplings hb 
and A 233 , respectively. These infrared fixed points provide a model independent theoretical 
upper bound on the B- violating coupling A 233 . It is worthwhile to point out here that the 
fixed point value for A 233 ~ 1.08 is more than a factor of 3 lower than what would be a naive 
perturbative upper bound of A 233 ~ ~ 3.5 for this coupling. 

In Figs. 4, 5 and 6, we present the fixed point behaviour of the corresponding A parame- 
ters. We notice the remarkable focussing property seen in the fixed point behaviour of all the 
A parameters. Again, we notice that the numerical evolution of the fixed point approaches, 
but does not actually reproduce, the approximate analytical values of Eqs. (|128|) , Q129D , 
and ( |13U| ). Since the quasi- fixed point value for the A parameter is inversely proportional 
to the quasi-fixed point value of the Yukawa coupling, it provides a lower bound on the 
corresponding A parameter. 
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V. SUMMARY AND CONCLUSIONS 



We have carried out a comprehensive study of the renormalization group flow in the 
nonminimal supersymmetric standard model with all the third generation Yukawa couplings, 
and the trilinear couplings A and k, and with highest generation baryon and lepton number 
violation. We have shown that the simultaneous fixed point for the top- and bottom- Yukawa 
couplings, and the B- violating coupling A233, is the only fixed point that is stable in the infra- 
red region. However, the top-quark mass predicted by this fixed point is incompatible with 
measured value of the top mass. This fixed point, therefore, provides a process- independent 
lower bound on the baryon number violating coupling A233 at the electroweak scale. 

We have shown that all other possible fixed point solutions are either unphysical, or 
unstable, in the infra-red region. In particular there is no infrared fixed point with simulta- 
neous B and L violation. This result is analogous to that found in MSSM with B and L 
violation [ |TT|j . 

We have also carried out the renormalization group analysis of the corresponding trilinear 
soft supersymmetry breaking parameters. We have obtained the true fixed points for these 
parameters, which serve as upper bounds on these parameters. 

Since the true fixed points are not reached in practice at the electroweak scale, we have 
also obtained the quasi-fixed points of the Yukawa couplings and the trilinear parameters. 
The quasi-fixed point values for the Yukawa couplings are numerically very close to the 
values obtained previously by ignoring the r Yukawa coupling. Since the quasi-fixed points 
are reached for large initial values of the couplings at the GUT scale, these reflect on the 
assumption of perturbative unitarity, or the absence of Landau poles, of the correspond- 
ing couplings. These quasi-fixed points, therefore, provide an upper bound on the relevant 
Yukawa coupling, especially the baryon number violating coupling A 2 3 3 . From the true fixed 
point and the quasi-fixed point ananlysis we are able to constrain the baryon number vi- 
olating coupling 0.97 ~ A233 ~ 1-08 in a model independent manner. We would like to 
emphasize that this infrared fixed point upper bound on A233 is more than a factor of 3 
lower than the naive perturbative upper bound on this coupling, and represents a significant 
constraint on it. We have complemented the quasi-fixed point analysis of the Yukawa cou- 
plings with an analysis of the corresponding soft supersymmetry breaking trilinear couplings. 
We have shown that the A parameters for the top- and bottom-quark Yukawa couplings, 
and the baryon number violating couplings all show striking convergence properties. This 
strong focussing property is quite independent of the input parameters at the unification 
scale (or equivalently the pattern of supersymmetry breaking), and the A parameters are, 
therefore, fully determined in the quasi-fixed regime. In particular, we have constrained 
the A parameters to be 0.77 ~ A t /nig ~ 1, 0.78 ~ A^jm-g ~ 1, and Ay /nig ~ 1. These 
constraints are analogous to those found in MSSM for these parameters. Thus, the infrared 
fixed point behaviour of NMSSM with B and L violation is similar to that of MSSM, and 
may be universal for supersymmetric models based on SU{2)l x U{1)y gauge group. 
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Figure Captions 



Fig. 1. The quasi-fixed point of the renormalization group evolution of the top-quark 
Yukawa coupling h t as a function of the logarithm of the energy scale. We have taken the 
initial values of h t at the grand unified scale M G ~ 10 16 to be 4.0, 3.0, 2.0, and 1.0, and 
evolved them down to ln(fi/M G ) ~ —33. The initial values of other Yukawa couplings are 
h b = 1.0, h T = 0, A = k = 0.5 and A'^ = 2.0. 

Fig. 2. The renormalization group evolution of the bottom-quark Yukawa coupling hb 
as a function of the logarithm of the energy scale. The initial values of h b at Mq are 
4.0, 3.0, 2.0, and 1.0. The initial values of other Yukawa couplings are h t = 1.0, h T = 0, A = 
k = 0.5, A' 2 ' 33 = 2.0. 

Fig. 3. Th equasi-fixed point behaviour of the baryon number violating Yukawa cou- 
pling A233 as a function of the logarithm of the energy scale. The initial values are 
A233 = 4.0, 3.0, 2.0, 1.0. The initial values of other other parameters are are h t = hb = 
1.0, hr = 0, A = k = 0.5. 

Fig. 4. Renormalization group evolution of ratio A t /m g as a function of the logarithm of 
the energy scale for several different initial values at the grand unified scale Mq. The initial 
values at Mq are —4, —3, —2, —1, 0, 1, 2, 3, 4. The initial values for other parameters at Mq 
are h t = 4.0, h b = 1.0, h T = 0, A = k = 0.1, X'^ = 1.0, and A b /m- g = 2.0, A x ^Jm g = 3.0. 

Fig. 5. Renormalization group evolution of the ratio Ab/m g as a function of the logarithm 
of the energy scale. The initial values of the ratio are — 4, — 3, — 2, — 1, 0, 1, 2, 3, 4. Other 
parameters at the scale M G are h t = 1.0, hb = 4.0, h T — 0, A = k — 0.1, A233 = 1.0, and 
A t /m- g = 2.0, A K Jm g = 2.57. 

Fig. 6. Renormalization group evolution of the trilinear coupling Ax'^ 33 /m g . The initial 
values of the ratio are —4,-3,-2,-1,0,1,2,3,4. Other parameters at the scale M G are 
h t = h b = 1.0, hr = 0, A = k = 0.1, and A t /m~ g = 2.0 A b /m g = 2.0. 
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FIG. 2. 
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34 




•35 -30 -25 -20 -15 -10 -5 



ln(n/Ny 

FIG. 6. 



35 



